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ABSTRACT 



Aims. We determine the age and mass of the three best solar twin candidates in open cluster M67 through lithium evolutionary models. 
Methods. We computed a grid of evolutionary models with non-standard mixing at metallicity [Fe/H] = 0.01 with the Toulouse-Geneva 
evolution code for a range of stellar masses. We estimated the mass and age of 10 solar analogs belonging to the open cluster M67. We made a 
detailed study of the three solar twins of the sample, YPB637, YPB 1 194, and YPB 1787. 

Results. We obtained a very accurate estimation of the mass of our solar analogs in M67 by interpolating in the grid of evolutionary models. 
The three solar twins allowed us to estimate the age of the open cluster, which is 3.87+jj 5| Gyr, which is better constrained than former estimates. 
Conclusions. Our results show that the 3 solar twin candidates have one solar mass within the errors and that M67 has a solar age within the 
errors, validating its use as a solar proxy. M67 is an important cluster when searching for solar twins. 

Key words. Stars: fundamental parameters - Stars: abundances - Stars: evolution - Stars: interiors - Stars: solar-type 
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1. Introduction 

Lithium depletion of each star depends on the evolution of 
its rotational history and the changes in the convective en- 
velope due to internal mixing. In main-sequence solar-like 
stars, lithium is easily destroyed by nuclear burning in stel- 
lar interiors at temperatures above 2.4 x 10 6 ^T, and its sur- 
face abundance indicates the depth of mixing below their pho- 
tospheres. The Sun has long been thought to be highly Li- 
depleted by a factor of 1 compared to the other field G stars 
dLambert & Reddvl2004. However, tw o of the best solar twins 
found bv lMelendez & Ramfrezl d2007l) (namely, HIP 56948 and 
HIP 73815) are Li -poor. The amount of Li depletion in solar 
twins is sensitive to microscopic diffusion, and an extra-mixing 
process is required to explain the low Li abundances observed, 
indicating that they also share a similar mixing history with the 
Sun. 

M67 is one of the most studied open clusters. Its chem- 
ical composition, metallicity and ag e are similar to those of 
the Sun. Various chemical analyses dTautvaisiene et al. hood 



Randich et al . 20061 IPace et aDl2008h show that not only abun 



dances of Fe, but also element abundances of O, Na, Mg, Al, Si, 
Ca, Ti, Cr, and Ni, are very similar to their solar counterparts, 
as close as allowed by the precision of the measurements. The 
ages determi ned are similar to those of the Sun, between 3.5 
and 4.8 Gyr dYadav et al.ll2008l) . Another interesting aspect of 



open cluster M67 is its Li-depleted G stars. Lithium abundance 
study is a crucial p oint for unders t andin g the physics involved 
in stellar interiors. IPasquini et al.l d 19971) show there is a large 
spread in Li abundances among solar-type stars in M67. Almost 
40% of their sample shows Li depletion comparable to that of 
the Sun. In this context, M67 is a perfect target to search for 
solar analogs and solar twins. Indeed, many stars are similar to 
the Sun in terms of effective temperature and chemical compo- 
sition, but their Li abundance is much higher t han in our star. 



Furthe rmore, among the 90 MS stars studied bv lPasquini et al 



( 2008) in M67, 10 had similar effective temperature and Li con- 
tent as low as that of the Sun, making them excellent solar twin 
candidates. 



In do Nascimento et al. I d2009l) . we showed that it is possi- 



ble to precisely determine the mass of solar twins using evolu 
tionary models. We determined t he mass and the age of 5 so 



lar t wins taken from two sources: [Melendez & Ramfrezl (|2007) 



and lTakeda et alJ (120071) . We began to compute a solar model, 
helioseismically calibrated, to reproduce not only solar lumi- 
nosity and radius, but also Li depletion. To that end, we intro- 



feet from the ^-currents ( 


Zahnlll992: Vauclair & Theadoll2003 


Theado & Vauclaii 2003 


). This calibration was used to com- 
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pute models of different masses whose tracks in a r e ff-log N(Li) 
diagram passed through the observational point. In the range 
of effective temperatures studied (typically 5400-6000 K), Li 
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Fig. 2. Same as Figure Q] for YPB1 194 



destruction is very sensitive to mass, and we were able to pre- 
cisely determine the mass of each solar twin and its age. 

In the present work, we computed helioseismically cali- 
brated, evolutionary, solar-like models with microscopic dif- 
fu sion and rotation-induced m ixing in the radiative interior, as 
d2009h . We estimate the masses of the 



do Nascimento et al 



m 

10 solar analogs in M67 found by iPasquini et al.l (120081) . We 
estimate the age of the three most probable solar twins of the 
sample. The working sample is described in Sect. [2] Improved 
Li abundances, using spectral synthesis instead of calibration 
based on equivalent width (EW, hereafter), are presented in 



Sect. [3] In Sect. |U we describe the grid of evolutionary mod- 
els with the non-standard physics that we computed. In Sect. [5] 
we present our results and compare our estimation of the age 
of M67 with earlier estimate. In this section we also show that 
the constitutive physics of the models influence and alter the 
isochrones constructed. We present the M67 color-magnitude 
diagram constrained by main sequence solar twins. Finally, our 
conclusions are outlined in Sect. [6] 

2. Sample of solar analogs in M67 

Pasquini et al. ( 20081) observed 90 targets belonging to cluster 
M67 with the multi-object FLAMES/GIRAFFE spectrograph 
(see Sect. [3]). They computed the effective stellar temperatures 
of these stars using two spectroscopic methods: the line-depth 
ratios (LDR) and the Ha wings. They derive d Li abundances 



from m easured EW using the calibration of ISoderblom et al 



(119931) . The authors highlight 10 stars in their sample, for 



which both T^° R and are within 100 K of the solar val- 
ues and which have a strong Li depletion. These stars are good 
solar twin candidates. However, as observed by the authors, 
given the limited resolution (R~ 17000) and S/N (80-1 10/pixel 
in average), the errors in Li abundance determination are high 
and the stars with upper limits may have a Li abundance com- 
parable to the solar value, or even higher. 

For our study, we used the effective temperature determined 
by the LDR method , which is believed to be more precise 



(IPasquini et al 



licity of [Fe/H] 
an age of 4.0 



20081) . For cluster M67 we adopted a m etal- 
0.01 + 0.03 dRandich et al.l l2006h and 



0.5 Gyr ( Dinescu et alJ 1995 ). Pasquini et al 



(2008) give a distance modulus for the cluster of (m - M)q 



9.63 + 0.06 sta t ± 0.05 sys , which corresponds to a distance d = 
843^* pc. Intrinsic absolute magnitudes My were derived from 
this distance, the V magnitude given by IPasquini et al.l ( 120081) . 
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and a visual extinction A y ~ 3.1 * E(B — V) = 0.127 mag 
dSchultz & Wiemei 1 1975b with E (B - V) = 41 ± 4 mmag, the 
reddening of M67 dTavlorll2007h . We computed relative stellar 
luminosity co mpared t o sola r luminosity, using the bolometric 
corrections of Flowei ill 9961) and the associated error from the 
uncertainty in distance. 

The atmospheric parameters of all these stars are in excel- 
lent agreement w ith those of the Sun. However, the Li deter- 
mination given in iPasquini et al. (2008) could be improved by 



spectral synthesis (see Sect. 13.11 ). The parameters of the solar 
analogs in M67 are summarized in Table Q] 

3. Observations and lithium abundance 

The spectra were acquired on three observation nights with 
the multi-object FLAM ES/GIRAFFE spectr ograph at the 
UT2/Kueyen ESO-VLT jPasquini et all 12002b in MEDUSA 
mod J] as part of a p roject searching for solar twins in M67 
(IPasquini et al 1 l2008h . The setting used was HR15N, which si- 
multaneously covers the Ha and Li i resonance doublet at 670.8 
nm with a resolution of R~17000. Three separate exposures 
were obtained to identify short and intermediate period bina- 
ries by comparing the radial velocities at different epochs. On 
average, combined spectra had a typical S/N of 80-1 10/pixel. 
For an extended descri ption of sample selec tion, observations, 
and data reduction, see IPasquini et al.l (12008b . 



models on the range of effective temperatures uncertainties. 
The procedure gives us an error of about 0.05 dex for log N(Li) 
due to uncertainty in effective temperature. The errors arising 
from uncertainties surface gravity and microturbulent velocity 
are less than 0.01 dex, so can be ignored. The total error in 
log N(Li) o btained by summing the error in model-fitting pro- 
cedure (see lBonifacio et al.l2007r) and effective temperature un- 
certainties in quadrature is equal to 0.10 dex as represented in 
Figs. H and 

4. Stellar evolutionary models 

For the purposes of this study, stellar evolutionary models 
were computed using the Toulouse-Geneva stellar evo- 



lution code TGEC dHui-Bon-Hoal 12008b . D etails on the 
physics of t hese models can be found in iRichard et al 



19961 12004. Ido Nascimento et alJ J2000t). and iHui-Bon-Hoc 



(2008). Standard input physics, non standard processes, 



diffusion, and rot ation-induced mixing added in the models 



described in do Nascimento et al. d2009l) . The method 



assumes the existence of an expected dependency between 
log jV( Li) and the star's mass and age f or solar analogs stars 
(e.g., Montalban & Schatzmanl 2000t Charbonnel & Talonl 



2005; 



Xiong&Dend 120091: Ido Nascimento et al.l I2069F 



3.1. Lithium abundance 



The Li abundance of all M67 twins was derived from the Li 1 
resonance transition at /16707 A synthetic spectrum fitted to 
the GIRAFFE spectrum for the fundamental atmospheric pa- 
rameters. First, we used T e ff spectroscopic determination and 
surface gravity \ogg = 4.44 dex, metallicity [Fe/H] = 0.0, 
projected surface velocity vsini = 2. km.s" 1 , mic r oturbu - 

lence velocity £ = 1.21 km.s 1 from IPasquini et al.l ( 12008b . 4.1. Models and calibration 



Baumann et alJ l2010b . This implies that stars with larger 
convection zones on the main sequence and a higher degree 
of differential rotation between the envelope and the radiative 
core present enhanced Li depletion, consequently we should 
know the angular momentum history of each studied star. 
Somehow, for solar twins stars belonging to the same cluster, 
their angular momentum average may be approximated well by 
a mean angular momentum evolution history. In this context, 
our analysis is limited to solar-like stars on the main sequence 
and with a similar angular momentum history. 



Atmospheric models were interpolated in the Kurucz grid 
jKuruczll 1993b. and the synthetic spectra were computed with 
MOOG dSnedenl [l973b . The GIRAFFE instrumental broaden- 
ing profile was set at 0.39 A, which is the FWHM of the 
GIRAFFE spectra. From the fundamental parameters, we de- 
rived logJV(Li) = 1.5, logJV(Li) < 1.3 and log JV(Li) = 1.6 re- 
spectively for YPB637, YPB1194, and YPB1787 (Fig. [Q Fig. 
12 and Fig.O, on the usual scale where logA^(H) = 12.00. 

Given the limited S/N and resolution of the YPB1194 
spectrum, only an upper limit for the Li abundance of this 
star can be determined. The consequences are discussed in 
Sect. |5l These value s are higher than the determinations of 
Pasquini etail d2008b . logjVYLi) = 1.4 and 1.0, and lower than 
0.6. Our Li abundance in the 10 solar twin candidates dif- 



We computed a grid of hundreds of evolutionary models of 
masses in the 0.90 to 1.10 M Q range with a step of 0.001 M , 
and with metallicities in the -0.030 to 0.050 range with a step of 
0.010, from the zero-age main sequence (ZAMS) to the end of 
the hydrogen exhaustion in the core. Our evolutionary models 
were calibrated to match the observed solar effective tempera- 
ture, luminosity, and Li abundance at the solar age. The cal ibra- 
tion method of models is based on the Richard et al. Idl996b pre- 
scription and is described in do Nascimento et al. "d2009b . In the 
following, we use models with a metallicity of [Fe/H] = 0.01, 
which is a simple average of d ifferent estimates of M67 metal- 
licity (see lPasquini et al.ll2008b . 



fered s omewhat from the estimates obtained by IPasquini et al.l 5. Results and discussion 



(2008) because we used spectral synthesis instead of a calibra- 
tion based on EWs, which may not be adequate for Li-depleted 
stars. For each star we determined logA^(Li) from synthetic 



1 This is the observing mode in FLAMES in which 132 fibers feed 
the GIRAFFE spectrograph. Some fibers are set on the target stars and 
others on the sky background. 



For each star we use an interpolation program to determine 
the mass. The uncertainties in T^ R give rise to uncertain- 
ties in mass determination. The mass is interpolated from an 
HR diagram and from a r e jf-log iV(Li) diagram. The results for 
each star are presented in Table[2] Determination from the T e g- 
log N{Li) diagram is more precise because Li abundance, more 
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Fig.4. Solar twin YP B1194 of open cluster M67 observed by 



Pasquini et al.l d2008l) : Li destruction along the evolutionary 
tracks as a function of effective temperature. The shaded zone 
represents the range of masses of TGEC models limited by the 
observational error bars. A model (dot-dashed lines) of 1.000 
M passes through the observed point. The position of the Sun 
and the lithium destruction of a solar model (dotted line) are 
indicated. The Li abundance of the solar twin is an upper limit 
(triangle). 



than luminosity, depends on mass (see Ido Nascimento et al 
20091) . All stars of the sample have a mass within the range 



expected for the mass of a solar twin (+5 % of the solar mass). 

In the following, we concentrate our analysis on objects 
YPB637, YPB1194, and YPB1787, which are the three most 
probable solar twins of the sample. They have a mass, an ef- 
fective temperature, and Li abundance very close to the solar 
values. The temperature errors in the Table[2]are those given in 
Table A2 of iPasquini et al.l (12008b . which are about 65K. 

For these three stars, we calculated an evolutionary model 
of the mass determined by interpolation in the r e ff-log N(Li) 
diagram, and we determined the range of masses limited by the 
observational error in T^ R , which gives a value of 1.005^ ^ 
M for the object YPB637, and 1.004+°°* M for the ob- 
ject YPB1787. For the object YPB1194, the determination of 
the Li abundance as an upper limit implies an upper limit for 
the mass determination: 1.000_o.oo8 M . In Figure|4]we present 
the r e ff-log A^(Li) diagram for YPB1194 as an example with 
the evolutionary tracks. The tracks of the evolution of the mod- 
els pass through the observation point, as expected. The posi- 
tion of the observation point on the track and the tracks limited 
by the error give us an age of 4.42+[ ^ Gyr for YPB637 and 
3.21+| ^ Gyr for YPB1787. For YPB1194, the upper limit of 
the Li abundance gives a lower limit for the age: 3.87 +2 03 Gyr. 
In Figure[5]we present Li destruction in the evolutionary model 
as a function of the age. 
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Fig. 5. Solar twins YPB6 37, YPB1194 and YP B1787 of open 
cluster M67 observed bv lPasquini et al.l d2008l) : Li destruction 
along the evolutionary tracks as a function of age. The positions 
of the solar twins, the Sun, and Li destruction of a solar model 
(dotted line) are indicated. The Li abundance of the solar twin 
YPB1 194 is an upper limit (triangle). 



5.1. Age determination of open cluster M67 

The usual way to determine the age of an open cluster is to 
fit a theoretical isochrone with the observed main sequence 
and turn-off stars in a color-magnitude diagram. The position 
of the observed points in the diagram depends on the adopted 
interstellar reddening for the cluster, and the position of the 
isochrone varies according to metallicity and distance modu- 
lus. The "hook" of the turn-off is re produced by calibrat i ng the 
overshooting of the convective core. iHobbs & Thorburnldl99lh 
determined an age of 5.2 + 1.0 Gyr from the high-dispersion 
echelle spectra of five stars l ocated along the main sequence of 
M67. In Montgomery et al.l d 1993k . the authors present a CCD 
photometry survey of M67. They derived a reddening E(B - 
V) = 0.05±0.01 and a metallicity [Fe/H] = -0.05+0.03. From 
cluster fitting to theoretical isochrones in a color-magnitude di- 
agram, they found an age between 3 a nd 5 Gyr, but none of 
the isochrones fit t he data consistently . Dinescu et al.l d 19951) 
used the sample of Montgomery et al.l d 19931) . cleaned up by 
cr oss-identifica t ion wi th the proper motion membership study 
of iGirard et alJ (Il989l) . They determined a distance modulus 
and the age of open cluster M67 by fitting the observed main 
sequence for single stars to the theoretical evolutionary tracks. 
They estimate an age of 4.0 + 0.5 Gyr, and found that earlier 
models with oversh oot of at most 0. 1 H p seem to reproduce the 
observations better. I Yadav et al.ld2008l) find an age between 3.5 
and 4.8 Gyr, by comparing their CMD, established from two- 
epoch archival observations with the Wide-Field Imager at the 
2.2 m MPG/ESO telescope and reduced with new astrometric 
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Fig. 6. Color-magnitude diagram of the solar analog sample of M67 from lYadav et al] (120081) . The dot-dashed isochrone, the 
continuous isochrone and the dashed isochrone correspond to an age of 3.21 Gyr, 3.87 Gyr, and 4.42 Gyr, respectively. The three 
best solar twins, YPB637, YPB1194, a nd YPB1787, are represented respectively as square, diamond, and triangle. The other 
solar analogs from lPasquini et al] ( 12008b are represented with open circles. 



techniques, as described in Anderson et al. ( 20061) . to the fol- 
lowing s ets of theoretical i sochro nes: the BaSTI models pro- 
vided by IPietrinferni et al] d2004 with E(B - V ) = 0.02 and 
[FeJH] = 0.06, the Padova models published by iGirardi et al. 
hood) with E(B -V) = 0.02 and [Fe/ H] = 0.00, the Yonsei- 



Yale isochrones provided by lYi et al] (1200 II) in their second 



release with E(B - V) = 0.02 and [Fe/H] = 0.00, and the 
Victoria-Regina models published bv lVandenBerg et al. I d2006l) 
with E(B - V) = 0.03 and [Fe/H] = -0.03. They confirm 
that the convec tive core overshoo t in M67 should be small, as 
pointe d out bv ISand auis tl d2004h and IVandenBerg & Stetson 



j2004 . IVandenBerg et al](l2007l) investigated the possibility of 



using the morphology of the M67 turn-off to put constraints on 
solar metallicity. Their models and isochrones computed on the 
assumption of a low-meta l mix b ased on the solar abundances 
derived by lAsplund et al.l (120051) . and adopting a metallicity 
[Fe/H] = 0.00, a reddening E(B -V) = 0.038, and a distance 
modulus (m - M)o = 9.70, do not predict the turnoff gap in the 
C MP. The same analy s is usin g a 3.9 Gyr isochrone assuming 
the Grevesse & Sauvall ( 1998 ) mix of heavy meta ls provides a 



good m atch to the morphology of the M67 CMD. Ma gic et al 



( 2010h confirm these results, but found that other physics in the 
models, e.g. element diffusion, nuclear reactions, prescription 
of core overshooting, also influence the stellar mass at which 



convective cores start to develop, and alter this result to the ex- 
tent that isochrones constructed with models using low CNO 
solar abundances can also reproduce the turn-off morphology 
in M67. They find ages between 3.9 and 4.8 Gyr. 

For each solar twin, we determined an age from the com- 
puted stellar model with mass previously determined. We 
found three different ages in a narrow range of 1.2 Gyr. The 
age that best fits the color-magnitude diagram is 3.87 Gyr, the 
upper limit of the age of YPB1194, the best solar twin can- 
didate (see Sect. 15.2b . We used the two other age determina- 
tions to determine an error bar. Hence, our estimation for the 
age of open cluster M67 is 3.87*2 gf Gyr, which is in very 



good agreement with the most recent results (IVandenBerg et al 
20071: lYadav etal]|2008l:lMagic et al]|2010b . 



5.2. The color-magnitude diagram 

Figure |6]presents the color-magnitude diagram of open cluster 
M67. The isochrones are constructed from computed models of 
metallicity [Fe/H] = 0.01 and masses from 0.90 to 1.31 M . 
The models of mass 1.19 M Q and higher are computed with a 
convective core overshooting parameter a ov = 0.01 Hp, whose 
occurrence is need ed to provide a sati sfactory match to the ob- 
served CMDs fsee lKalirai et alj|200U and references therein). 
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A very small convective core overshoot is needed to match 
the feature s near the main s e quence turnoff region, in agree - 
me nt with ISandauistl (|2004 . IVandenBerg & Stetson! d2004 . 
and lYadav et al.l (12008b. We used equation (3) and Table 4 of 
Casagrande et al.l d2010l) to compute the (B - V) color index 
of our models. We determined the absolute visual magnitude 
of the models from the luminosit y using (Mboi)p - 4.74 mag, 
and the bolometric corrections of lFlowerl (11996b . For the stars 
of the M67 sample, w e used the (B - V) and V measurements 



of lYadav et al.l (120081) . Our best fit of the main sequence and 
the "hook" of the turn-off is achieved with isochrone 3.87 Gyr, 
using a distance modulus (m - M)q = 9.68 and a reddening 
E(B - V) — 0.02 mag. The value of the dist ance modulus is in 
good agreement with former determinations dVandenBerg et al 



2007tlYadav et al.l2008l:|Pasquini et al.l2008b. T he reddening is 
lower t han the valu e obtai ned bv lTavlorl d2007l) but, as pointed 
out by lYadav et al.l (I2008fl . given the intrinsic uncertainties in 
calibration of E(B - V), the true reddening could be anywhere 
between 0.02 and 0.04 mag. 

6. Conclusions 



In this paper, we appl ied the method developed by 
do Nascimento et al. d2009l) to the sol a r anal ogs of open clus- 



ter M67 identified by iPasquini et a l. (2008), in order to de 



termine their masses as accurately as possible. We identified 
three objects, YPB637, YPB1 194, and YPB1787, whose effec- 
tive temperature, luminosity, mass, and Li abundance values 
are very close to the solar ones. For these three solar twins 
we derived ages of 4.42+}-*j Gyr for YPB637, and 3.21+[ ^ 
Gyr for YPB1787. For YPB1194, the lower limit of the age 
is 3.87 +203 Gyr. These estimations allowed us to construct 
isochrones on a color -magnitude d i agram for open cluster M67, 



using the sample of lYadav et al.l (|2_008). The isochrone 3.87 



Gyr matches both the main-sequence and the turn-off regions, 
using a distance modulus (m - M)q = 9.68 and a reddening 
E(B - V) = 0.02 mag. The three solar twins allowed us to esti- 
mate the age of the open cluster as 3.87^° ^ Gyr. We thus con- 
firm the importance of M67 as a cluster with a slightly lower 
age than the solar one, and the mass of the solar twin candi- 
dates is also very similar to solar; therefore, M67 is an impor- 
tant cluster when searching for solar twins. 
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Table 1. Parameters of the observed solar analogues in M67 . 
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T LDR 
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(K) 


(L/L Q ) 


\ogN(U) NLlh 


logAT(Li) moo / 
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1 1 12 +01!9 

1 - llz '-0.106 


<0.6 
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1 089 +0 117 
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1.5 
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5756 ± 60 


r, q 7 c+0.099 

U.7Z.J_ 008g 


<0.6 


< 0.6 


1194 


5766 ± 64 


976 +0 105 


<0.6 


< 1.3 
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5716 + 64 


0.9603$ 


<0.6 


1.2 


1304 


5704 + 64 


0.886!°;^ 


<0.5 


< 0.8 


1315 
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1.8 
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5716 + 63 


n ST 1 +0.088 
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1,UJ -0.099 


1.0 


1.6 
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<0.5 
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1 Note: Cols 1-4 from Pasquini et al. (2008). 
Our determinations of the Li abundances using the MOOG program 



Table 2. Results of mass determination from interpolation in the r e ff-log N(hi) diagram. 
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